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Abstract 

Polarized Compton scattering on the deuteron is studied in nuclear effective field theory. A 
' set of tensor structures is introduced to define 12 independent Compton amplitudes. The scalar 

Q ■ and vector amplitudes are calculated up to 0{{Q/ A)'^) in low-energy power counting. Significant 

CN ■ contribution to the vector amplitudes is found to come from the spin-orbit type of relativistic 

^ . corrections. A double-helicity dependent cross section Aia = ((t+i„i — cr+i+i) /2 is calculated to 

the same order, and the effect of the nucleon isoscalar spin-dependent polarizabilities is found to 
be smaller than the effect of isoscalar spin-independent ones. Contributions of spin-independent 
polarizabilities are investigated in various asymmetries, one of which has as large as 12 (26) percent 
^ ■ effect at the center-of-mass photon energy 30 (50) MeV. 
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I. INTRODUCTION 



Compton scattering is an important tool to probe the internal structure of a composite 
system, such as an atomic nuclei. As quantum electrodynamics involved in the process 
is well understood, the remaining uncertainty is associated with the strong interactions 
among nucleons in nuclei. Thus Compton scattering data enables physicists to extract 
information about the nuclear structure from the underlying strong interaction dynamics. 
Recent progress in high-energy, high-intensity photon beams has made Compton scattering a 
practical tool for nuclear physicists [l| . In particular, a polarized photon beam is capable of 
studying spin aspects of strong interaction physics. This paper focuses on polarized Compton 
scattering on the deuteron, the double-helicity dependent cross section in particular, in the 
framework of nuclear effective field theory (EFT). 

The deuteron, as the simplest nuclear system, is of great importance for understanding the 
nucleon-nucleon interactions and the properties of individual nucleons. Polarized Compton 
scattering on the deuteron presents a new opportunity to probe spin physics. Indeed, because 
the deuteron is a loosely bound system, one might expect to learn a host of spin- dependent 
properties of the neutron and proton as free particles. This possibility is especially important 
for the structure of the neutron because there is no free neutron target in nature. 

It has been realized for some time that nuclear physics at low energy might be under- 
stood by effective field theories (EFT) which works according to the same principles as the 
standard model 12]. However, constructing a workable scheme for specific systems is not 
straightforward. In the past few years, considerable progress has been made in the two- 
nucleon sector (see |3] for a recent review). It began with the pioneering work of Weinberg, 
who proposed to encode the short-distance physics in a derivative expansion of local oper- 
ators [2]. The problem associated with the unusually small binding energy of the deuteron 
was solved by Kaplan, Savage and Wise by exploiting the freedom of choosing a renor- 
malization substraction scheme 0], quickly followed by the pionless version ^ (see also 
[1; IHl)- Required reproducing the residue of the deuteron pole at next-to- leading order 
(NLO), a version with accelerated convergence was suggested in p|. The use of dibaryon 
fields as auxiliary fields, first introduced in was taken seriously in 11 1 which simplified 
the calculation significantly. 

From the viewpoint of nuclear EFT, Compton scattering on the deuteron at low energy 
can be divided into two regions according to the photon energy u. Region I is where the 
photon energy is far below the binding energy of deuteron B = 2.2 MeV and, hence, u/B 
is a small parameter. Region II is where the photon energy is above the binding energy, 
but significantly below the mass of the pion, for example, ~ 50 MeV. In Region I, one 
makes the low energy expansion of Compton araplitudes and studies various polarizabilities 
of the deuteron defined through the expansion [12]. Studies in this ultra- low energy region, 
where the binding effect plays a dominant role, provide insight about the internal structure 
of the deuteron as a bound state. In Region II, the probing photon is more sensitive to the 
responses from individual nucleons. Therefore, Compton scattering there may serve as an 
alternative tool to study free-nucleon properties, such as spin-independent and dependent 
polarizabilities. In this paper, we are mostly interested in the second region. 

Extracting the isoscalar spin-independent polarizabilities ao and /3o from unpolarized 
Compton scattering has attracted considerable attention in the past two decades. Although 
there are three types of amplitudes (scalar, vector and tensor) contributing to the cross 
section, only the scalar amplitudes have been included in some of the calculations of the 
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unpolarized cross section so far. Nuclear EFT seems to provide a justification for this. 
However, a recent work 3| showed that vector amphtudes contribute significantly (order 
of 15% or more) to the unpolarized cross section, because of the enhancement from a factor 
of the square of the isovector magnetic moment /x^. It turns out that this enhancement 
has its effect not only on unpolarized scattering but also on polarized one, leading to, 
for instance, a bigger helicity-dependent cross section. Although this makes it easier to 
measure it experimentally, the effect also diminishes the contribution from the isoscalar 
nucleon spin-dependent polarizabilities, and hence makes it harder to access them from the 
future Compton data. 

To demonstrate the above point, we calculate a double- helicity dependent (vector- 
polarized) cross section up to the order at which the spin polarizabilities contribute, and 
compare the results with and without their contribution. The photon-nucleon interactions 
considered in this calculation include the electric current and magnetic couplings, and the 
spin-orbit terms from the non-relativistic reduction of the relativistic interactions. It has 
been realized previousl y th at the relativistic corrections are surprisingly large in potential 
model calculations EE [ll . In EFT, the spin-orbit interactions were taken into 

account in the studies of the deuteron forward spin-dependent polarizabilities and Drell- 
Hearn-Gerasimov sum rule jl^. They were neglected in other EFT calculations because they 
are nominally suppressed in power counting by l/M^ relative to the other two couplings. 
However, for certain spin-dependent observables, their contributions can be of leading order, 
as we shall see. 

The paper is organized as follows. Section II is devoted to kinematics, where we write 
down 12 basis structures for scattering amplitudes using parity and time-reversal symmetries. 
The scalar and vector structures are the same as those in Compton scattering on a spin-1/2 
particle such as the proton. The tensor structures are new, and useful for general discussions 
of polarized deuteron Compton scattering. Section HI explains a calculation of the vector 
Compton amplitudes using the dibaryon formulation of EFT. Power counting in both Re- 
gions I and II is explained to show the significant contribution of the spin-orbit interactions. 
The result of individual diagrams is listed in Appendix B. In Sec. IV, a double-helicity- 
dependent (vector-polarized) cross section is defined, and the numerical result is presented 
with and without the contribution from the nucleon spin-dependent polarizabilities. The 
feasibility of using polarized Compton data to extract these polarizabilities is discussed. In 
Sec. V, we investigate the effect of the spin-independent polarizabilities on a number of spin 
asymmetries. Section VI contains the conclusion of the paper. 



II. REAL PHOTON-DEUTERON COMPTON SCATTERING AMPLITUDES 

In this section, the general tensor structure of the amplitudes for real photon Compton 
scattering on a deuteron is considered. Through helicity counting, it is easy to see that there 
are a total of 12 independent amplitudes. We choose these amplitudes in a basis convenient 
for subsequent calculations. We comment on the frame dependence of the tensor structures 
associated with the amplitudes. 

The real photon has two independent helicities ±1; the deuteron has three, ±1 and 0. 
Therefore, the total number of helicity amplitudes is 2x3x2x3 = 36. Parity invari- 
ance of strong and electromagnetic interactions restricts the number of independent ones 
to 36/2=18. Among those, time-reversal symmetry relates 6 to the others with initial and 
final state exchanged. This reduces the number of independent amplitudes to 18 — 6 = 12. 
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Moreover, the general result of helicity counting can be derived, and is 2(J+ 1)(2J+ 1) for 
a spin-J target. 

In the low-energy region, it is convenient to use the nonrelativistic notation for tensor 
structures associated with the amplitudes. If the spins of the initial and final deuterons are 
coupled, the sum is 0, 1 or 2. The amplitudes classified in this way are called scalar, vector, 
and tensor, respectively. Clearly the number of scalar amplitudes must be the same as that 
of Compton scattering amplitudes on a spin-0 target, namely, 2; and the number of vector 
amphtudes is the same as that on a spin-1/2 target, 4. Thus the number of independent 
tensor amplitudes is 12 — 2 — 4 = 6. 

In the remainder of this section, we construct a set of 12 linearly-independent structures, 
using the 3-momenta of the photon and deuteron, and their polarization vectors. Among 
four 3-momenta, only three are independent because of the momentum conservation. By 
choosing a specific frame, one more constraint follows, and hence only initial and final 

— * — * 

3-momenta of photon, k and k', are needed for the construction. The initial and final three- 
momenta of the deuteron, p and ff can be expressed in terms of these of the photon. For 
example, the lab frame is defined by p = and = k — k', the center-of-mass frame (CM) 
hj p = —k and p' = —k', and the so-called Breit frame by p = — k) and = —-^{k' — k) 
and so p + = 0. The constraints among momenta associated with a frame are generally 
not invariant under symmetries such as time reversal, which exchanges the initial and final 
momenta and reverses their directions, and photon crossing symmetry, which exchanges the 
initial and final photon with the sign of energy and 3-momentum flipped. For instance, the 
momentum constraint in the lab frame is not invariant under either time reversal or crossing 
symmetry, while the momentum constraint in the CM frame violates crossing symmetry. 

According to the above, in the CM and Brcit frames where parity and time reversal 
invariance are manifest, there are 12 independent tensor structures for the Compton ampli- 
tudes. These structures are constructed out of initial and final photon polarization vectors 
(e'* and e), deuteron polarization vectors (^'* and ^), and the initial and final photon mo- 
mentum vectors k — k/\k\ and k' — k' /\k'\. One can couple (^'* and ^) into scalar, vector 
and tensor to obtain scalar, vector, and tensor amplitudes. Alternatively, these structures 
can be obtained by the matrix element of a unit matrix /, spin matrices Jj, or tensor 
( JjJj -|- JjJi — trace) between the initial and final deuteron polarization vectors. Under par- 
ity transformation, all momentum and polarization vectors change sign, whereas the spin 
matrices do not. Under time-reversal transformation, these quantities transform according 
to: e ^ e'*; k <^ -k'; J ^ -J. 

Requiring symmetry under both parity and time-reversal, we choose the 12 basis struc- 
tures for Compton scattering on the deuteron as follows. 

Pi = e'* ■ e / = e'* ■ e e'* ■ e , 

P3 = iJ.e'*xe = (rxO-(e'*xe), 

P4 = i J . s'* X S = X ■ (s'* X S) , 

P5 = i{J ■ y* -e-J- k'e'* ■ s) = [{C X ■ ks'* ■ e - x ■ k'e'* ■ s] , 
P6 = i{J- k's'* -e-J- he'* ■ s) = [{t x ■ k's'* ■ e - {C x ■ ke'* ■ s] , 

^ _( J . u ■e'* + J- e'*J ■ e - ^e'* ■ e) = t ■ e'*i ■ e + t ■ ei ■ e'* - '^t ■ ie'* ■ e , 
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P8 = -{J ■ sj ■ s'* + J ■ s'*j ■ s - ^s'* ■ s) = i'* ■ s'*i ■ s + i'* ■ si ■ s'* - '^i'* ■ is'* ■ S , 

pg = -e'* . k{J .eJ-k' + J-k'J-e)-e- k'{J ■ e'* J ■ k + J ■ kJ ■ e'*) + -e'* ■ ke ■ k' 

3 

= e'* ■ Hi'* ■ k'i ■ e + i'* ■ ei ■ k') + e ■ k'{i'* ■ ki ■ e'* + i'* ■ e'*i ■ k) - ^i'* ■ ie'* ■ ke ■ k' , 
pio = -s'* .k(j.sJ-k' + J-k'J-s)-s- k'{J ■ s'*J -k + J-kJ ■ s'*) + \s'* ■ ks ■ k' 

= s'* ■ k{i'* ■ k'i ■ s + i'* ■ si ■k') + s- k'ii'* ■ ki ■ s'* + i'* ■ s'*i ■ k) - y,'* ■ ir -ks-k' , 

pn = -e* -eiJ-kJ-k + J- k'J ' k' ~ ^) = e* ■ e{i'* ■ ki ■ k + i'* ■ k'i ■ k' - " , 

Pi2 = -r •s(j-fcj-^ + j-^v-^'-^) = r -^'e-^'-p'-o , (i) 

where the s and s'* are defined as s = A; x e and s'* = k' x e'*. These structures are constructed 
in such a way that duahty between the electric and magnetic fields is manifest. Under the 
dual transformation, e ^ s, s ^ —e, which is a 7r/2-rotation in the photon polarization, 
the above structures transform as p2i-i P2i with i = 1,...,6. The structures with the 



unit matrix and spin operators (pi to pe) are the same as those for a spin-1/2 target j21 
Appendix A explains why these 12 structures are complete and independent. 

The most general Compton scattering amplitude on the deuteron can be expressed as 

/ = E f^P^ , (2) 
1=1 

where fi defines the spin-dependent amplitudes. The first two {i = 1,2) are scalar ampli- 
tudes; the following four {i = 3, ...,6) are vector amplitudes; and the last six {i = 7, 12) 
are tensor amplitudes. 



III. VECTOR COMPTON AMPLITUDES TO (Q/A)'^ FROM EFT 

In this section, we calculate the vector Compton amplitudes to (Q/A)'^ in a low-energy 
expansion in nuclear EFT. The calculation is based on the dibaryon approach in the pionless 
theory, which has been referred as dEFT{ji) 11]. 

A central concept in EFT is power counting. EFT is designed to describe physics at 
one scale — low-energy scale in this case — using an effective lagrangian, and the physics at 
other scales is accounted for through the couplings. Power counting allows a systematic way 
to take into account corrections from other energy scales. For Compton scattering on the 
deuteron, the natural momentum scale is ^/M^B (Mat is the nucleon mass) which will be 
generically referred to as Q. The deuteron binding energy B is then counted as order of Q^. 
The energy and momentum of the external photon probe, = is counted as 

• m Region I where u <^ B, and as 

• Q in Region II, where u ~ y/Mj^B. 
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The high-energy scales include the nucleon mass Mjv, the pion mass mjr, and similar scales 
describing the structure of the nucleon, like the charge radius, and parameters in nucleon- 
nucleon interactions. Because and are very different, we use A to denote scales at 
around m^, and identify rriT^/Mi^f as Q/A. Therefore, ratio Q/Mj^r can actually be treated 
as (Q/A)"^. Although this is not fully consistent in the EFT sense, it is a way to organize 
numerically-close ratios phenomenologically j2ol |. 

In dEFT(/r), the nucleon rescattering in both singlet ^5*0 and triplet ^5*1 channels is 
represented by the propagation of dibaryon fields, tj and Sa, respectively. The lagrangian 
density for the triplet channel is 
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where A^ is the 2-component nucleon field with an implicit isospin index. The time and 
spatial derivatives with electromagnetic gauge symmetry are Dq and D, respectively. = 
-^r20"2crj is the projection operator of the triplet channel, and y is the coupling between 
nucleons in the triplet channel and the triplet dibaryon. Requiring producing the nucleon- 
nucleon scattering amplitude, one has 



y 



in 



A 



(4) 



with ij, being the renormalization-scale introduced in the power divergent subtraction scheme 
0]. The parameters a and r are the scattering length and effective range, respectively. In the 
present formulation, these two are counted as order in both singlet and triplet channels. 
Thus the scaling property of y and A is ?/ ~ \/Q and A ~ Q^, respectively. Dressing the 
dibaryon propagator with nucleon bubbles does not change the counting of the propagator. 
Therefore the bubbles must be summed to all orders; the dibaryon propagator dressed with 
nucleon bubbles is 



An 



H + 



47r 



(A-E 



+ i\ MnE 



(5) 



with E the center-of-mass energy. The wave function renormalization constant is the residue 



at pole E = —B, and a simple calculation yields [ll|: Zd = 7r*^ ~ ^•') . 

We remark that it is straightforward to convert the nuclear EFT lagrangian with the 



nucleon field into that in dEFT(7f). Following the prescription in [ll|, one converts a pair 



of nucleon fields in the singlet and triplet channels to dibaryon fields. 



N^P'^N 



(6) 



where = :^o"2T2Ta is the projection operator for the singlet channel. 

Nuclear EFT describes the interactions between the nucleons and external electromag- 
netic probes systematically. Besides the coupling generated in the covariant derivatives in 
the above lagrangian density, D = V — zeA, there is also the magnetic coupling to the 
nucleon. 



2M, 
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-Art(/io + /iir3)(T-BAr 



(7) 
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where no and /ii are the nucleon's isoscalar and isovector magnetic moments, respectively. 
An associated term is the spin-orbit-type relativistic correction 



^so = ((2/io + (2^1 - I) ra) iVt^ • (D x E - E x D)iV , (8) 

which is generated from the reduction of a relevant relativistic interaction. 

There are also interaction terms involving the dibaryon fields themselves. One term 
accounts for the transition between the ^5*1 and ^5*0 channels through a magnetic field, 

-^^em 1 = e , / , =t\s-iBj + h.c. . (9) 

The coupling constant Li has been determined by the rate of n + p — (i + 7. The measured 
cross section a = 334.2±0.5 mb with an incident neutron speed of 2200 m/s fixes Li = —4.42 
fm. Another term involves the elastic scattering of the deuteron in the magnetic field, 

£em,2 = (/io - ^) s'^'iB^t, , (10) 

with the value of L2 fixed to be —0.03 fm from the magnetic moment of the deuteron. The 
/in in the above equation is introduced to reproduce the magnetic moment at leading order 
22 . 23! . There is also an associated relativistic correction. 



^-'^ =2k:(^'~A~^) '"''^ (D X E - E X D)^. t, , (11) 

which generates a seagull interaction of the dibaryon and electromagnetic fields. At last, 
there are nucleon polarizabilities interactions: 

£poi = 27tN^ («o + am) E^N + 27tN^ {(3o + /3iTs) B^N 

+27cN^ (7^ + 7^x3) CT-EXEN + 2nN^ (7^ + 7^X3) a-BxBN 
-2nN^ (7^ + 7^x3) EijCTiBjN + 2nN^ (7^ + 7S^3) B,,a,E,N , (12) 

where Eij = l/2(ViEj + VjEj) and Bij = l/2(VjBj + VjBj) are the electric and magnetic 

field gradients. The nucleon isoscalar {(yo,l3o,'JEi,Mi,E2,M'2) ^''^^ isovector (ai,/9i,7£^i,Mi,£;2,M2) 
polarizabilities are defined as, for example, ao = 1/2 (op -|- a„) and ai = 1/2 (ap — an), with 
similar relations for others. The isoscalar ones are what can be probed in deuteron Compton 
scattering. Chiral perturbation theory calculations yield 23 |: 

ao = lOpo = 12 X lO^W, 

7^ = -3.1 X 10-Hm\ 7S = 0.4 x 10-^fm^ 

7^ = 2.1x10"^, 7^?^ = 0.6 X lO^W. (13) 

Feynman diagrams that contribute to the deuteron Compton scattering to [Q/A)^ in 
power counting are shown in Figs. 1-4. Figure 1 contains diagrams with direct photon- 
dibaryon interactions. Figure 2 contains the seagull interactions with the nucleon. The 
diagram 2c actually corresponds to the contribution from electromagnetic polarizabilities 
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(a) - Q' 



FIG. 1: Compton scattering with photons coupled to the dibaryon field directly. The open circle 
denotes the electric photon- dibaryon coupling from the gauged derivative. The solid dot denotes the 
seagull term in Eq. The intermediate thick line represents the triplet dibaryon. 






FIG. 2: Diagrams with seagull interactions on the nucleon lines. The small open circle denotes the 
coupling from the gauged derivative in the first term in Eq. the small solid circle represents 
the coupling from spin-orbit interaction defined in Eq. while the small open box represents the 
point interactions associated with polarizabilities of the nucleon in Eq. Power counting of 

the leading contribution of each diagram is listed below the diagram. In diagram 2c, the two count- 
ings, Q and Q^, are for spin-independent and spin- dependent nucleon polarizabilities contributions, 
respectively. 



of the nucleon. Figure 3 include diagrams without intermediate dibaryon fields. Finally 
diagrams in Fig. 4 have intermediate singlet and triplet dibaryon propagations. 

To estimate the importance of a particular diagram in our power-counting scheme, we 
need to study the dominant regions of a loop momentum in the integral. Let us use 
q) to denote the loop momentum generically. The size of the loop momentum is de- 
termined by poles of the propagators. Typical nucleon propagators in the loop integration 
are i/ [—B — Qo — 2m '^hen the photon momentum does not pass through the nucleon 



4jy -(t^H 4:^ 

(e)~Q2 if)-Q'/' {g)-Q"' [h)^Q^I^ 

FIG. 3: Diagrams without intermediate dibaryons. The small open circles denote the electric 
photon-nucleon coupling from the gauged derivative in the first term in Eq. the small shaded 
circles denote the magnetic photon-nucleon coupling in Eq. while the small solid circles repre- 
sent the spin-orbit interaction between photon and nucleon in Eq. 
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FIG. 4: Diagrams with intermediate dibaryon states. The small open circles denote the electric 
coupling in Eq. and the small shaded circles denote the magnetic coupling in Eq. The 
intermediate thick lines with one arrow represent both the spin singlet and triplet channels. The 
solid box denotes the Li and L2 couplings in Eqs. ^ and il i^) . 

line, and i/ (^qo + 00 — ^ + ie^ when the photon momentum does. Because go scales as 



the former has a momentum pole at |g| ~ vi? and the latter a pole at |g| ~ y/u = y\k\. In 
Region I, these two poles have the same order of magnitude and have power counting |g| ~ Q. 
In Region II, the pole (|(?| ~ VB) has counting \q\ ^ Q, and the other pole ~ has 
|g| ~ ^/Q. A Feynman integral can be approximated by the pole that produces a leading 
contribution. 

For example, let us count the power of diagram (b) in Fig. 3. The Feynman integral 
has a momentum power Quj'^\q\^/{\(^^uj), where Q is from the wave function renormaliza- 
tion, tu^ is from two magnetic couplings, the u in the denominator is from the propagator 
V (q'o + ^ ~ ^ + ) ^iid |g| is the loop momentum, with d'^q counted as and three 
other propagators in the denominator as At the pole |g| ~ Q, it is of order Q; and 
at the other pole \q\ ~ y/Q, it is Q^^'^. Thus the leading contribution is of order Q, shown 
below the diagram. Note that the Q-counting here is dimensionally balanced by the nucleon 
mass Mtv in the denominator. 

Because there are multiple leading regions in a Feynman diagram, power counting can be 
rather tricky sometimes. For example, the nominally higher-order, spin-orbit couplings 
can produce leading contributions in a certain momentum region. To see this, let us 
compare the power counting for diagrams (f) and (h) in Fig. 3. The counting of (f) is 
Q\q\^{q + k)uj / {\q\^uj'^) , where in the denominator \q\'^ is from the two propagators that do 
not depend on the photon momentum and u;^ is from two propagators that do; in the nu- 
merator {q + k) and u factors are from the derivative and magnetic couplings, respectively. 
Since only the k term in {q + k) survives the symmetrical momentum integration, diagram 
(f) is of order Q^^"^. On the other hand, counting of diagram (h) is Q\q\^{q + fc)^co'/(|g|^a;^) 
which, compared to diagram (f), has an extra power of {k + q) / Mj^, because it is a relativistic 
correction. However, the dominant term contributing to the integral is (f in the {q + kY 
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factor, which is of order Q at the leading pole. Therefore, diagram (h) is also of order Q^/^. 
Thus the spin-orbit coupling contributes as significantly as the magnetic coupling in these 
diagrams. 

Power counting allows us to determine the leading contribution of every Feynman di- 
agram. The result is indicated below each diagram in Figs. 1, 2, 3, and 4. Again, all 
countings so far are in terms of powers of Q/M^^ including that for the nucleon polariz- 
ability in diagram 2c. We will treat them as (Q/A)^ in phenomenology as we mentioned in 
the beginning of this section. According to chiral perturbation theory, the spin-independent 
polarizabilties contribute to the scalar amplitudes at order (Q/A)^ j20|, the spin- dependent 
ones contribute to the vector amplitudes at order (Q/A)'^. An explanation of the count- 
ing of the nucleon polarizability contributions from diagram 2c is in order. Compared 
with the leading order contribution T2a in Appendix B, the result of T2c is suppressed by 
2Mm{q.ooj'^ 1 Pouj'^, 7^] jyfi E2 M2^^) /<^em, which is numerically (Q/A)^ for scalar polarizabilties 
and (Q/A)^ for vector ones. 

According to the above, the scalar amplitudes start at (Q/A)^, vector amplitudes at 
(Q/A)^ and tensor amplitudes at {Q/A)^. However, the leading-order vector amplitudes are 
actually proportional to the square of the isovector magnetic moment nf, and are numerically 
larger than what simple power counting indicates. Therefore, the vector-vector contribution 
to the unpolarized cross section is quite significant On the other hand, the enhancement 
makes the contribution of nucleon spin-dependent polarizabilties relatively less important. 

From Figs. 1 to 4, the vector-polarized amplitudes can be calculated to order (Q/A)^. Our 
explicit results are shown in Appendix B. In order to have the result look more compact, 
the integration over Feymann parameter x has not been completed. One must exercise 
caution, however, when the power of the un-integrated result is counted. For example, the 
result of diagram (b) in Fig. 3 seems to scale as 'ju'^ / (Mj^icu)^^'^ ~ Q^^'^- However, after the 
integration, it actually scales as Q, consistent with power counting. 

IV. A DOUBLE-HELICITY DEPENDENT (VECTOR-POLARIZED) CROSS 
SECTION 

With the scalar and vector amplitudes presented in the previous section and Appendix 
B, we can calculate spin-dependent Compton scattering cross sections. Of course, any 
polarized cross section can be constructed out of the complete 12 (scalar, vector, and tensor) 
amplitudes once they are known. Because the tensor amplitudes start at order (Q/A)^, we 
do not need to know them to predict certain spin-dependent cross sections up to some orders 
in Q. 

As we have seen in the previous section, the vector amplitudes receive contribution from 
the spin-dependent polarizabilities of the nucleon. Therefore, we would like to find a cross 
section which can be used to probe the vector amplitudes, and hence possibly extract the 
spin polarizabilities. 

A double-helicity dependent cross section satisfies the above condition. Suppose the 
helicities of the initial-state photon and deuteron are Ai and Ai, respectively. The general 
Compton scattering cross section with these polarized initial states is o"aa- Define a vector- 
polarized cross section 

^if^ = 2 (^+1-1 ~ ' (14) 

where -|-1 (—1) is a right-handed (left-handed) polarization. If the initial momentum of the 
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FIG. 5: The vector-polarized cross sections for different CM frame photon energy u) = 30, 50, 70, 90 
MeV. (See text for comments on the 70 and 90 MeV cases). 6 is the scattering angle in the 
CM frame. The dashed lines contain no contribution from spin-independent or spin- dependent 
polarizabilities. The dotted lines have contributions from spin-independent polarizabilities of the 
nucleon, but without dependent ones of nucleon. The solid lines have contribution from both. The 
values of nucleon polarizabilities are taken from chiral perturbation theory in Eq. 



photon is along the z direction, the scattered photon momentum is taken along a direction 
with a polar angle 6. Then the polarization vector of the in-coming photon is e = — ^(ic + 
iy). The deuteron, moving in the negative z direction with a negative helicity, has the same 
wave function. The deuteron with a positive helicity has a wave function ^ = — iy). 
Note that the beam is circularly-polarized in the so defined vector-polarized cross section. 
Actually, investigations indicate that the vector amplitudes cannot be probed as leading 
order contributions if the beam is parallel-polarized. 

According to the above definition, the vector-polarized Compton cross section can be 
expressed in terms of the full 12 amplitudes as follows, 

Aia = Re[S*V + V*V + V*T + T*T] (15) 



2M2 



N 



3{u + Ju^ + MlY 



-Re 



[-6(1 + z^){flf, + /*/4) - i2z(A74 + fih 
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-6^(3 + z'){rj, + f;fe) - 6(1 + Sz'Xrje + f^h)] 

+[-3(1 - z'){f;fs + flh + nu + flh) 

-zzii - z^)u;h + flh)] + [-4(2 - z^){f;h + flh) - Mf^h + flfr) 

-5z{l - z'){f;h + flfm) - 5(1 - z^){fU,o + flf9) + (1 + 7z^)(/37ii + flfi^) 
+z{h + ?>z^)UU\2 + flfii) - ^(9 - z')U;f-j + /*/8) - (5 + ?>z'){f;h + /e/r) 
-2(1 - z%hh + /eVio) - 4^(1 - ^')(/57io + /eVg) + 2z(3 + 5^^)(/57ii + /eVis) 
+ (1 + 12^2 + 3/)(/57i2 + /e/ii)] + [-3(3 + z^){f;fj + /^/s) + ^Azf^h 
+^z{l - z')U;h + /s/io) - 15(1 - z^)U;f\o + /s/g) + 3(1 - z^)U; fii + f;fi2) 
+3^(1 - z^){f;fn + /gVii) - 6(1 - z^)\f;h + /i*o/io) 

+6;^(1 - z^){f;fn + /ro/12) + 3(1 - z%f;f,2 + /ro/11)]] , (16) 

where 5"*^, V*V , V*T, and T*T denote combinations of scalar- vector, vector- vector, vector- 
tensor, and tensor-tensor amplitudes, respectively. According to power counting, the domi- 
nant contribution is from the scalar and vector interference, and is of order (Q/A)^. If cal- 
culating the cross section to order (Q/A)^, we need the scalar amplitudes to order (Q/A)^ 
and vector amplitudes to (Q/A)^, including the nucleon polarizability term. The tensor 
amplitudes do not contribute at this order. Therefore, the vector-polarized cross section is 
a useful observable to probe the vector amplitudes, and hence the spin polarizabilities. 

We have shown in Fig. 5 the vector-polarized cross section to {Q/KY in EFT at CM 
photon energy u — 30, 50, 70, 90 MeV, respectively. The contribution from spin-independent 
polarizabilities of the nucleon is more significant at higher energy. There is virtually no 
difference between the cross sections with the polarizabilities turned on or off at the photon 
energy a; = 30 MeV. However, there is a notable difference at 50 MeV and substantial 
difference at 70 MeV and 90 MeV. [Note, however, that our results for 70 and 90 MeV 
are just for exploratory study, because the pion has to be included as a dynamical degree 
of freedom at such high energies. However, we expect that the general features will not 
change in a full analysis.] The effect of the nucleon polarizabilities is more significant at 
forward and backward angles (almost zero at 9 = 7r/2). Moreover, the contribution from 
spin-independent polarizabilities ao,/3o is of similar size at forward and backward angles, 
while the spin- dependent polarizabilities contribute mainly at forward angles. 

According to power counting, both the scalar and spin polarizabilities contribute to the 
vector-polarized cross section at order (Q/A)^. However, the leading-order vector amplitude 
is enhanced by a factor nf. Therefore, the scalar polarizabilities contribute more significantly 
to the cross section, and generate a larger infiuence than the spin polarizabihties. As seen in 
the figure, AiO" — especially at the backward angles — is very sensitive, as is the unpolarized 
cross section, to the scalar polarizabilities of the nucleon. Therefore one cannot extract the 
vector polarizabilities without knowing the scalar ones to a reasonable accuracy. From Fig. 
5, the best way to extract the spin polarizabilities is to measure Aicr at forward angles and 
at relatively high energy (higher than 50 MeV). On the other hand, the EFT expansion 
becomes less reliable at high energy. 
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V. ASYMMETRIES SENSITIVE TO SPIN-INDEPENDENT NUCLEON POLAR- 
IZABILITIES 



As seen from the previous section, the spin-independent nucleon polarizabihties have 
to be determined before the extraction of spin- dependent ones become possible. In this 
section, we investigate various asymmetries with the goal of extracting spin-independent 
polarizabihties. 

Asymmetries are generally easier to measure than cross sections because of the cancella- 
tion of systematic errors. The asymmetry associated with the vector-polarized cross section 
in the previous section is: 

= (17) 

where the indices ±1 have the same meaning as in Eq. ()14|) . The expression for the 
numerator has been shown in the previous section. The expression for the denominator in 
terms of scalar and vector amplitudes is: 

1 , 

2 l^+i-i + 



2M 



-Re 



;i + ^')(/i7i + /272) + 4^(/i72 + nf^) + 2(/373 + rj^) 



+1(3 + 12^2 ^ ^4)(^*^^ ^ ^ 2;.(5 + 3z2)/576 + (3 + ^z^mu + nu) 

+z{7 + z'){f;f, + f:U)] (18) 

The result of for CM photon energy u = 30, 50 MeV is shown in Fig. 6. Clearly, 
as the vector-polarized cross section, the asymmetry at the backward angle has stronger 
dependence on oq, compared to other angles and shows almost no sensitivity on spin- 
dependent polarizabihties. However, unlike the cross section, the dependence on the ao, 
in the asymmetry is suppressed to about 8% at 50 MeV due to cancellation between the 
numerator and denominator. 

In the following, we investigate other asymmetries in search of a larger dependence on 
ao, Po- There are two new asymmetries related to when the polarization axis of the 
deuteron target is changed. If the xz plane is chosen as the scattering plane, one can define 
an asymmetry with deuteron polarized linearly in the x direction: 

with the first index +1 of a indicating that the photon is right-handed polarized, the second 
index indicating that the deuteron target is polarized in the Jx = ±1 states. The expressions 
for the numerator and the denominator in terms of scalar and vector amplitudes are: 

2 A/f 2 

^ VT^Re \-2 (zflh + flh + (1 + z^)flh + 2^A76 + fUz 



+zf;f, + 2zf;h + (1 + z^)f;h) + zu;h + flh) + 2/374 + (1 + z')f;h + 2^/376 
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180 



FIG. 6: The asymmetry for different CM frame photon energy io = 30, 50 MeV. 9 is the 
scattering angle in the CM frame. The meaning of dashed lines, dotted lines, and solid lines are 
the same as in Fig. 5. The values of nucleon polarizabilities are taken from chiral perturbation 
theory in Eq. 111,^) . 
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FIG. 7: The asymmetry T,^ for different CM frame photon energy u) = 30, 50 MeV. 6 is the 
scattering angle in the CM frame. The meaning of dashed lines, dotted lines, and solid lines are 
the same as in Fig. 5. The values of nucleon polarizabilities are taken from chiral perturbation 
theory in Eq. \lt!\) . 



+2zflf, + {l + z')flf, 



2Ml 



-Re 



;i + z'){fUi + /2/2) + 4^/r/2 + (2 - z'){f;f, + fif,) 



+2zf:j, + (3 + z'){f;f, + f*j,) + z{5 - z'){f;f, + fu,) 



(20) 



The result for S^. at CM photon energy uj = 30, 50 MeV is shown in Fig. 7. The peak of 
this asymmetry is around the scattering angle of 105 degrees, where the dependence on a^, 
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FIG. 8: The asymmetry for different CM frame photon energy uj = 30, 50 MeV. 9 is the 
scattering angle in the CM frame. The meaning of dashed lines, dotted lines, and solid lines are 
the same as in Fig. 5. The values of nucleon polarizabilities are taken from chiral perturbation 
theory in Eq. 111,^) . 

Po is about 8% at 50 MeV. 

Similarly, one can define the asymmetry with the deuteron polarized in the y direction, 
which is perpendicular to the scattering plane. It turns out this asymmetry is actually a 
single-spin asymmetry, independent on the polarization of the photon beam. 

= ^"-=+^ ~ ^"-=-^ , (21) 

where the photon beam is unpolarized and the deuteron target is polarized in the Jy = ±1 
states. The expressions for the numerator and the denominator in terms of scalar and vector 
amplitudes are: 



^ K=+i - ^j.=-i) = - \ '^^]^^y^^ ^i^[fih + f;f3 + zf*js + zf;f,] 



2 ^ ' ' ^ (uj + Juj'^ + Mj.}^ 



(i + z'){f*f, + f;f,) 

+1(2 - z') u;h + fih) + + 3^') uih + f;h) 

+(1 + z^) u;h + flh) + '2z (A72 + nh + flh) 

+zf;f, + z{3 + z')f;f, (22) 

The result for Hy at CM photon energy = 30, 50 MeV is shown in Fig. 8. The peak of 
this asymmetry is around a scattering angle of 90 degree, where the dependence on oq, /3o 
is about 12% at 30 MeV and 26% at 50 MeV, much larger than the dependence in E^^. 
Therefore, the single-spin asymmetry should serve as a good observable to extract nucleon 
scalar-isoscalar polarizabilities. Note that the polarizations of the deuteron in the above 
asymmetries are defined in the CM frame, while in experiment the deuteron is prepared 
polarized in the lab frame. The polarization in these two frames are different in case of ^x,y 
This is an error of size uj/Md which can be safely neglected at low energy. 
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The tensor amplitudes contributions are not taken into account in these results shown 
above. They are small contributions from the analysis of the power counting. But numeri- 
cally, the effect of them could be enhanced due to the large size of isovector nucleon magnetic 
moment, which also explains that the vector amplitudes effect are enhanced. While a more 
complete calculation of asymmetries with all the tensor amplitudes included is beyond the 
scope of this paper, we did, however, study their effects on the as ym metries by using the 
tensor amplitude ff from a previous calculation in EFT with pion j25|. We found that the 
Sj^ is less dependent on these amplitudes compared with the other asymmetries, which offers 
an additional reason that this asymmetry is better than others for extracting ao and jSo- 

We have also investigated the parallel-perpendicular single spin asymmetry, which is the 
ratio of the difference and sum of two cross sections when the deuteron target is unpolarized 
and the photon beam is linearly polarized either parallel or perpendicular to the scattering 
plane. This asymmetry is found to have a weak dependence (about 3% at 50 MeV) on qq, 
Po than and therefore is not presented here. 

VI. CONCLUSION 

In this paper, we presented a convenient set of basis for Compton scattering on the 
deuteron. We then calculated the scalar and vector Compton amplitudes to C((Q/A)'^) 
in a nuclear EFT without the pion, at which the scalar and spin polarizabilities of the 
nucleon contribute. The result was then used to calculate a double-helicity-dependent cross 
section which is linearly proportional to the vector amplitudes. We studied the effects of 
the polarizabilities on the cross section, finding that the scalar polarizabilities have more 
dominant influence than the spin polarizabilities. Thus an accurate measurement of the 
cross section can help to determine the former. However, if the scalar polarizabilities are 
determined with good accuracy, the cross section can provide a constraint on the spin- 
dependent ones. Finally, we investigated various asymmetries in search of large dependence 
on scalar polarizabilities and found that T^y has the best potential. 

This work was supported by the U. S. Department of Energy via grant DE-FG02-93ER- 
40762 and by the National Science Council of Taiwan, ROC. JWC thanks Paulo Bedaque 
for organizing the Summer Lattice Workshop 2004 at Lawrence Berkeley Laboratory where 
part of this research was completed. 

APPENDIX A: TENSOR BASIS FOR DEUTERON COMPTON AMPLITUDES 

The 12 basis structures can be systematically obtained by keeping track of the matrix 
structure sandwiched between the initial and final deuteron polarization states. The struc- 
tures with unit matrix and single spin matrix are the same as the structures for spin-1/2 
target. There are six such structures (pi ~ pe) [HI- Our goal is to find out the remaining 
six structures, which should all be of tensor type with symmetrized double spin matrixes. 

To write them down, first notice that since there are double Js associated with them, the 
parity invariance requires that there are an even number of cross products among vectors: 
e,e'*,k,k' and two Js. Moreover, since any even number of cross products can be trans- 
formed into dot products, we only need to write down structures with dot products. Since 
subtracting trace is straightforward, we choose to do it at the end. The structures before 
subtracting trace can be found systematically by looking at which pair dot with Js and what 
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is left over. First, if the pair is e and e'*, There is only one such structure: 

Ti = J-eJ-e'* + J-e'*J-e . (Al) 
If the pair is k and k', there are two structures: 

T2 = e'* ■ e{J ■ kJ ■ k' + J ■ k'J ■ k) , 

Tg = L e'*k' -eiJ-kJ-k' + J- k'J ■ k) . (A2) 

If the pair is k and e'*, time reversal invariance requires that the other pair k' and e appear 
in the same structure and in the following combination: 

T4 = e'* ■k{J-eJ-k' + J-k'J-e)+e- k'{J ■ e'*J ■ k + J ■ kJ ■ e'*) . (A3) 

If the pair is k and e, time reversal invariance rcqiiircs that the other pair, k' and i'*, appear 
in the same structure and in the following combination: 

Ts = e'* -kiJ -eJ -1^ + J -kJ ■e) + e- k'{J ■ e'*J -k' + J-k'J- e'*) . (A4) 

If the pair is two ks, the time reversal invariance requires that the other pair, two k's, appears 
in the same structure and in the proper combination. There are two structures of this type: 

Tg = e'* -eiJ -kJ ■k + J ■ k'J ■ k') , 

Tr ^ k- e'*k' -eiJ-kJ-k + J-k'J- k') . (A5) 

The above way of constructing structures with double Js exhausts all possibilities. There 
is no problem about the completeness. However, we get more structures than expected from 
helicity counting. It is hard to find the relation among them directly and it turns out that 
we need to make use of the duality character of the electric magnetic field. Starting from 
the above seven structures, we can write down another set of structures which covers the 
above set and has the duality correspondence among them, just like the structures from pi 
to pe- Without knowing the dependence among the structures from ti to ry, the minimal 
number of such a set of structures is eight. They are chosen as: 



= J-eJ-e'* + J-e'*J-e, 
= J-sJ-s'* + J-s'*J-s , 
= e'* .e(J-kJ-k' + J- k'J ■ k) , 
= s'* ■ s{J -kJ -k' + J ■ k'J ■ k) , 

= e'* .k(J.eJ-k' + J-k'J-e) + e- k'(J ■ e'*J ■ k + J ■ kJ ■ e'*) , 
= s'* -kiJ-sJ-k' + J -k'J ■§) + §■ k'(J ■ s'*J -k + J-kJ- s'*) , 
r; = e'* -eiJ-kJ-k + J- k'J ■ k') , 

§'* ■ s(J -kJ ■k + J ■ k'J ■ k') . (A6) 



One notices that under duality transformation, these eight structures transform as: T2j_i -v^ 
T2j, with i — 1,2, 3, 4. This set with eight structures can be expressed in term of seven /jS 
and the expression is found to be: 

t[=Ti , T2= 4p2 + T4 - ZTi - Ta , T3 = , T4 = ZT2 - T3 , 

T5 = T4 , Te = 2zT2 + T5- 2t3 - 2tq , Tj ^ , Tg = ZTe - T-j . (A7) 
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Since eight structures are expressed in terms of other seven structures, one relation among 
r/s must exist, and it is found to be: 



zr[ + r'^ + r'^-r'^ = Ap^ , (A8) 

from which another relation can be found through duality transformation of the above 
relation: 

zr'^ + r[ + T'^-T'^ = 4pi . (A9) 
Now, we have two constraints on eight structures and are therefore left with six independent 
structures, as expected from helicity counting. We choose t( 2,5,6,7,8 basis structures. 

With trace subtracted explicitly, they are: 

P7 = -(J . ej . e'* + J . e'V ■ e - ^e'* • e) = i'* ■ e'*i ■ e + i'* ■ ■ e'* - '^i'* ■ ^e'* ■ e , 

P8 = -{J ■sJ-s'* + J- s'*J ■ s - U'* ■ s) = r • s*i ■ s + i'* ■ si ■ s'* - '^i'* ■ is'* ■ s , 

P9 = -e'* ■ k{J -iJ -k' + J -k'J -e) -e- k'{J ■ e'* J ■ k + J ■ kJ ■ e'*) + ^e'* ■ ke ■ k' 

= e'* ■ k{i'* ■ k'i ■ e + i'* ■ ei ■ k') + e ■ k'{i'* ■ ki ■ e'* + i'* ■ e'*i ■ k) - -i'* ■ ie'* ■ ke ■ k' , 

3 

= -s'* .k{J -sJ -k' + J -k'J ■s)-s- k'{J ■ s'*J -k + J-kJ- s'*) + -s'* -ks-k' 

3 

= s'* ■ k{i'* ■ k'i ■ s + i'* ■ si ■k') + s- k'ii'* ■ ki ■ §'* + i'* ■ s'*i ■ k) - ^i'* ■ is'* -ks-k' , 

Pn = -e* .e{J-kJ-k + J- k'J ■k'~^) = e'* ■ e^'* ■ H ■ k + i'* ■ k'i ■ k' - ^i'* ■ i) , 

= -s'* ■ s{j ■ kj ■ k + J ■ k' J ■ k' - \) = s'* ■ s{i'* ■ki-k + i'* ■ k'i ■ k' - ■ i) -(Aio) 

PiS = 1 ~ 12) are the basis structures of deuteron Compton amplitudes in the frame 
where time reversal invariance is manifest such as the Breit frame and center-of-mass frame. 
Note that lab frame is not such a frame because it lacks the symmetry between the initial 
and final deuteron. 

There are other tensor structures that are often met in studies of Compton scattering on 
the deuteron. Here we provide a list and their relation to the basis set defined above: 

e'* ■ e{i'* ■ ki ■ k' + i'* ■ k'i ■ k) = -zp, - P8 + P9 + \zpi , 

e'* ■ ke ■ k'ii'* ■ ki ■ k' + i'* ■ k'i ■k) = il- z')pr + zp, - pio + ^z'pi - \zp2 , 

2 

e* ■ ke ■ k'{C ■ k^ ■ k + ^ ■ k'^ ■ k') = zpu - pu + -{zpi - P2) , 

e'* . k(i'* .ei.k + i'* ■ ki ■ e) + e ■ k'{i'* ■ e'*i ■ k' + i'* ■ k'i ■ e'*) = 2p, + 2zp^ - pio + 2pn , 

e'*-e(rxO-(A:'x^) = 2;p3 + p4-p5, 

(r X i) . (e'* X k)e ■ k' - ii'* X i) . (e x k')e'* ■ k = 2zp, - p, , 

ii'* X i) ■ se'* ■ k - ii'* X i) . s'*e ■ k' = 2p, - p, , (All) 
where z = cos 9 = k -k', which is used throughout this paper. The last three expressions for 



the vector- type structures have appeared in the literature before [2l|. Other useful relations 



can be obtained from the above through the duality transformation. 
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APPENDIX B: COMPTON AMPLITUDES TO (Q/A)^ IN EFT 

Diagrams with the photon directly coupled to the dibaryon are shown in Fig. 1. The 
result is: 
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a;p3(l-4/.o + ^) • 



(Bl) 



Diagrams with the seagull interaction on the nucleon line are shown in Fig. 2, among 
which are contributions from nucleon polarizabilities. The result for each diagram is: 
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(B2) 



with associated with the nucleon polarizabilities. 

The contribution without the intermediate singlet or triplet state is from diagrams in Fig. 
3. The result of each diagram along with photon crossing and the diagram with interchange 
of two photon coupling vertices, if different, is: 
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--(2pi + pe) / cia;-- )- ^ . ,,,, + (2pi - pe) dx- ^ ^ 



2V ' '^°'io (72 - M^^u;x - ie)3/2 ' '^"^ Jo {j^ + M^ux - ie)^/^ 

2 



X ( / dx . ^ - [ dx 
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+ Mr^cux — ie Jo Vl'^ ~ M^ujx — ie ) ' 



g2 



y dx\J 7^ + Mj^ujx — ie — J dx\J 7^ — Mj^ujx — ie^ . (B3) 

The diagrams with the intermediate triplet or singlet state are shown in Fig. 4. The 
result from each diagram along with photon crossing and the diagram with interchange of 
two photon coupling vertices, if different, is: 
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